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Abstract 

PsJ ■ In this paper, a multi-source multi-relay cooperative wireless network with binary modulation and binary 

_ network coding is studied. The system model encompasses: i) a demodulate-and-forward protocol at the 

. relays, where the received packets are forwarded regardless of their reliability; and ii) a maximum-likelihood 

optimum demodulator at the destination, which accounts for possible demodulations errors at the relays. An 

■ asymptotically-tight and closed-form expression of the end-to-end error probability is derived, which clearly 
_ showcases diversity order and coding gain of each source. Unlike other papers available in the literature, the 

,— — i, proposed framework has three main distinguishable features: i) it is useful for general network topologies 

■ and arbitrary binary encoding vectors; ii) it shows how network code and two-hop forwarding protocol affect 
^ ■ diversity order and coding gain; and ii) it accounts for realistic fading channels and demodulation errors at 

\ the relays. The framework provides three main conclusions: i) each source achieves a diversity order equal 

, to the separation vector of the network code; ii) the coding gain of each source decreases with the number 

0^ . of mixed packets at the relays; and iii) if the destination cannot take into account demodulation errors at 

■ the relays, it loses approximately half of the diversity order. 
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I ! 
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sJ ■ I. Introduction 

Cooperative communications and Network Coding (NC) have recently emerged as strong candidate tech- 
nologies for many future wireless applications, such as relay-aided cellular networks [1], [2]. Since their 
inception in [3] and [4], they have been extensively studied to improve performance and throughput of 
wireless networks, respectively. In particular, theory and experiments have shown that they can be extremely 
useful for wireless networks with disruptive channel and connectivity conditions [5]-[7]. 

However, similar to many other technologies, multi-hop/cooperative communications and NC are not 
without limitations [1], [8]. Due to practical hardware limitations, e.g., the half-duplex constraint, relay 
transmissions consume extra bandwidth, which implies that using cooperative diversity typically results in 
a loss of system throughput [9]. On the other hand, NC is very susceptible to transmission errors caused by 
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noise, fading, and interference. In fact, the algebraic operations performed at the network nodes introduce 
some packet dependencies in a way that the injection of even a single erroneous packet has the potential to 
corrupt every packet received at the destination [10], [1 1]. Due to their complementary merits and limitations, 
it seems very natural to synergically exploit cooperation and NC to take advantage of their key benefits while 
overcoming their limitations. For example, NC can be an effective enabler to recover the throughput loss 
experienced by multi-hop/cooperative networking, while the redundancy inherently provided by cooperation 
might significantly help to alleviate the error propagation problem that arises when mixing the packets [1]. 

In this context, multi-source multi-relay networks, which exploit cooperation and NC for performance 
and throughput improvement, are receiving an always increasing interest for their inherent flexibility to 
achieving excellent performance and diversity/multiplexing tradeoffs [12]-[36]. More specifically, consider- 
able attention is currently devoted to understanding the achievable performance of such networks when both 
cooperation and NC operations are pushed down to the physical layer, and their joint design and optimization 
are closely tied to conventional physical layer functionalities, such as modulation, channel coding, and 
receiver design [37], [38]. In particular, how to tackle the error propagation problem to guaranteeing a given 
quality-of-service requirement, e.g., a distributed diversity order, plays a crucial role when these networks 
are deployed in error-prone enviroimients, e.g., in a wireless context. For example, simple case studies 
in [16], [22], [23], and [39] have shown that a diversity loss occurs if cooperative protocols or detection 
algorithms are not adequately designed. To counteract this issue, many solutions have been proposed in 
the literature, which can be divided in two main categories: i) adaptive (or dynamic) solutions, e.g., [20], 
[23], [27], [28], and [35], which avoid unnecessary error propagation that can be caused when encoding 
and forwarding erroneous data packets; and ii) non-adaptive solutions e.g., [16], [19], [21], [22], and [25], 
which allow erroneous packets to propagate through the network but exploit optimal detection mechanisms 
at the destination to counteract the error propagation. Each category has its own merits and hmitations. 

Adaptive solutions rely, in general, on the following assumptions [23], [24], [27], [28]: a) network code 
and cooperative strategy are adapted to the channel conditions and to the outcome/reUabiUty of the detection 
process at the relay nodes. This requires some overhead, since the network code must be communicated to 
the destination for correct detection; b) powerful enough channel codes at the physical later are assumed to 
guaranteeing that the error performance is dominated by outage events (according to the Shannon definition 
of outage capacity) [36, Sec. II]; and c) the adoption of ideal Cyclic Redundancy Check (CRC) mechanisms 
for error detection, which guarantees that a packet is either dropped or injected into the network without 
errors (i.e., erasure channel model). However, recent results have shown that, in addition to be highly spectral 
inefficient as an entire packet is blocked if just one bit is in error, relaying based on CRC might not be 
very effective in block-fading channels [40], [41]. An interesting hnk-adaptive solution, which does not 
require CRC for error detection and avoids fuU-CSI (Channel State Information) information at the relays, 
has been proposed in [18]. Therein, the achievable diversity (using the Singleton bound) is studied under 
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the assumption that ad hoc interleavers are used, while no analysis of the coding gain is conducted. 

Non-adaptive solutions rely, in general, on the following assumptions [16], [19], [25]: a) neither error 
correction nor error detection mechanisms are needed at the physical layer, but the relays just regenerate 
the incoming packets and forward them to the final destination (i.e., error chaimel model). This results in 
a simple design of the relay nodes, as well as in a spectral efficient transmission scheme as the received 
packets are never blocked; and b) the possibiUty to receive packets with errors needs powerful detection 
mechanisms at the destination, which require CSI of the whole network to counteract the error propagation 
problem and to achieve full-diversity. Similar to adaptive solutions, this requires some overhead. 

As far as adaptive solutions are concerned, [23], [27], [28] have recently provided a comprehensive study 
of the diversity /multiplexing tradeoff for general multi-source multi-relay networks, and have shown that 
the design of diversity-achieving network codes is equivalent to the design of systematic Maximum Distance 
Separable (MDS) codes for erasure channels. Thus, well-established and general methods for the design of 
network codes exist, which can be borrowed from classical coding theory. On the other hand, as far as non- 
adaptive solutions are concerned, theoretical analysis and guidelines for system optimization are available 
only for specific network topologies and network codes. To the best of the authors knowledge, a general 
framework for performance analysis and code design over fading channels is still missing. Motivated by these 
considerations, in this paper we focus our attention on non-adaptive solutions with a threefold objective: 
i) to develop a general analytical framework to compute the Average Bit Error Probability (ABEP) of 
multi-source multi-relay cooperative networks with arbitrary binary encoding vectors and realistic channel 
conditions over all the wireless links; ii) to provide guidelines for network code design to achieve a given 
diversity and coding gain tradeoff; and iii) to understand the impact of the error propagation problem and 
the role played by CSI at the destination on the achievable diversity order and coding gain. 

More specifically, by carefully looking at recent literature related to the performance analysis and code 
design for non-adaptive solutions, the following contributions are worth being mentioned: i) in [16], the 
authors study a simple three-node network without NC (a simple repetition code is considered) and they show 
that instantaneous CSI is needed at the destination to achieve full-diversity. No closed-form expression of 
the coding gain is given; ii) in [19] and [33], the authors introduce and study Complex Field Network Coding 
(CFNC), which does not rely on Galois field operations and exploit interference and multi-user detection to 
increase throughput and diversity. The analysis is valid for arbitrary network topologies. However, only the 
diversity order is computed analytically, while the coding gain is studied by simulation; iii) in [21], the authors 
study a simple three-node network with binary NC. UnUke other papers, chaimel coding is considered in the 
analysis. However, the error performance is mainly estimated through Monte Carlo simulations; iv) in [22], 
the author considers multiple relay nodes but a simple repetition code is used (no NC). Main contribution 
is the study of the impact of chaimel estimation errors on the achievable diversity; v) in [25], the authors 
study a network topology with multiple sources but with just one relay. Also, a very specific network code 



3 



is analyzed. This paper provides a simple and effective method to accurately computing the coding gain of 
error-prone cooperative networks with NC; vi) in [34], the authors analyze generic multi-source multi-relay 
networks with binary NC, but error-free source-to-relay links are considered, and the performance (coding 
gain) is computed by using Monte Carlo simulations; vii) in [39] and [46], we have studied the performance of 
network-coded cooperative networks with realistic source-to-relay wireless channels. However, the analysis 
is useful only for two-source two-relay networks and for a very specific binary network code; viii) in [42], a 
general framework to study the ABEP for arbitrary modulation schemes is provided, but a simple three-node 
network without NC is considered; and ix) in [43], the authors study a three-node network with a simple 
repetition code. Exact results are provided for coding gain and diversity order. Finally, in [44] and [45], 
NC with error-prone source-to-relay links is studied, but the analysis is applicable only to noisy channels, 
while channel fading and distributed diversity issues are not investigated. 

According to this up-to-date analysis of the state-of-the-art, it follows that no general framework for 
performance analysis and design of non-adaptive solutions exists in the literature, which is useful for 
generic network topologies, for arbitrary encoding vectors, and which provides an accurate characterization 
of diversity order and coding gain as a function of the CSI available at the destination. Motivated by 
these considerations, in this paper we focus our attention on a general multi-source multi-relay network 
with realistic and error-prone channels over all the wireless links. For analytical tractabiUty (and to keep 
the implementation complexity of relays at a low level [34], [47]), we consider a binary network code, a 
Binary Phase Shift Keying (BPSK) modulation, and the Demodulate-and-Forward (DemF) relay protocol. 
With these assumptions, the main contributions and outcomes of this paper are as follows: i) a Maximum- 
Likelihood (ML-) optimum demodulator is proposed, which allows the destination to exploit the distributed 
diversity inherently provided by cooperation and NC. The demodulator takes into account demodulation 
errors that might occur at the relay nodes, as well as forwarding and NC operations. It is shown that 
the demodulator resembles a Chase combiner [48] with hard-decision decoding at the physical layer; ii) a 
simple but accurate framework to compute the end-to-end ABEP of each source is proposed. The framework 
provides a closed-form expression of diversity order and coding gain, and it clearly highlights the impact 
of error propagation and NC on the end-to-end performance; iii) it is proved that each source node can 
achieve a diversity order that is equal to the separation vector [49], [50] of the network code. In particular, it 
is shown that the optimization of network codes is equivalent to the design of systematic linear block codes 
for fully-interleaved fading channels, and that Equal and Unequal Error Protection (EEP/UEP) properties 
are preserved [49]; and iv) the impact of CSI at the destination is studied, and it is shown that half of the 
diversity order is lost if the destination is unable to account for possible demodulation errors at the relays. 

The paper is organized as follows. In Section JIJ network topology and system model are introduced. In 
Section JIIJ the ML-optimum demodulator that accounts for demodulation errors at the relays is proposed. 
In Section |IVj a closed-form expression of the end-to-end ABEP is given. In Section |Vl diversity order and 
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coding gain are studied for arbitrary binary network codes and network topologies. In Section |Vll numerical 
results are presented to substantiate analysis and findings. Finally, Section IVIII concludes this paper. 

II. System Model 

We consider a generic multi-source multi-relay network with Ns sources (St for t = 1,2, Ns), 
Nji relays (Rq for g = 1,2, . . . , A'^), and, without loss of generality, a single destination D. We consider 
the baseline Time Division Multiple Access (TDMA) protocol, where each transmission takes place in 
a different time-slot, and multiple-access interference can be neglected [3]. We assume that direct links 
between sources and destination exist, and that the relays help the sources to deliver the information packets 
to the final destination. The cooperative protocol is composed of two main phases: i) the broadcasting phase; 
and ii) the relaying phase. During the first phase, the source St transmits the information packet intended to 
the destination in time-slot for t = 1, 2, . . . , Ns- These Ns packets are overheard by the Nr relays too, 
which store them in their buffers for further processing. This phase lasts Ns time-slots. During the second 
phase, the relay Rg forwards a linear combination, i.e., NC is applied [4], of some received packets to the 
destination in time-slot Tj^^+g for g = 1, 2, . . . , Nr. We consider a non-adaptive DemF relay protocol, which 
means that each relay demodulates the received packets, but perform NC and forward them regardless of their 
reliability. As a result, packets with erroneous bits can be injected into the network. However, these packets 
can be adequately used at the destination, by exploiting advanced detection and signal processing algorithms 
at the physical layer, to improve the system performance [1, pp. 18-20]. According to the working operation 
of the protocol, broadcasting and relaying phases last Ns + Nr time-slots. Since Ns information packets 
are transmitted by the sources, the protocol offers a fixed rate, TZ, that is equal to 7^ = Ns/{Ns + Nji). In 
this paper, we are interested in understanding how the operations, i.e., NC, performed at the relays affect the 
end-to-end performance for this given rate. Main objective is understanding the performance of cooperative 
networks with NC when physical layer terminologies are exploited to counteract the error propagation 
problem [37], and, more specifically, when demodulation and network decoding are jointly performed at 
the destination (i.e., cross-layer decoding). For analytical tractability and simplicity, we retain three main 
reasonable assumptions: i) uncoded transmissions with no channel coding are considered. Accordingly, 
there is no loss of generality in considering symbol-by-symbol transmission. Some preliminary results with 
channel coding are available in [51]; ii) BPSK modulation is assumed to keep the analytical complexity at 
a low level; and iii) binary NC at the relays is investigated. However, unlike many current papers in the 
literature, e.g., [25], [39], [46], and references therein, no assumption about the encoding vectors is made. 
These assumptions are widespread used in related literature e.g., [14], [16], [22], and the references therein. 

A. Broadcasting and Relaying Phases 

According to the assumptions above, the generic source St broadcasts, in time-slot Tt, a BPSK-modulated 
signal, xst, with average energy Em, i.e., xs^ = \/£'m (1 — 265^), where bs^ G {0, 1} is the bit emitted by 
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St- Then, the signals received at relays ii^ for g = 1, 2, ... , and destination D are: 

yStD = hstDXSt + nStD 

where hxY is the fading coefficient from node X to node Y, which is a circular symmetric complex Gaussian 
Random Variable (RV) with zero mean and variance per dimension (Rayleigh fadingj). Owing to the 

distributed nature of the network, independent but non-identically identically distributed (i.n.i.d.) fading is 
considered. In particular, let dxY be the distance between nodes X and Y, and a be the path-loss exponent, 
we have = d^y [52], [53]. Also, nxY is the complex Additive White Gaussian Noise (AWGN) at 
the input of node Y and related to the transmission from node X to node Y. The AWGN in different time 
slots is independent and identically distributed (i.i.d.) with zero mean and variance Nq/2 per dimension. 

Upon reception of ystR^ and ystD in time-slot Tt, the relay Rg for q = 1,2, . . . , Nr and the destination 
D demodulate these received signals by using the ML-optimum criterion, as follows: 



bStRq = argmin <! \ystRg - -/E^hs^R^ (l - 2bst 



bstO = argmin <j \ystD - VE^hs^D (1 - 26s, 



(2) 



arg mill 

6St6{0,l} 

where (•) denotes the demodulated bit and (•) denotes the trial bit used in the hypothesis-testing problem. 
More specifically, bstR^ and bs^D are the estimates of bs^ at relay i?^ for g = 1, 2, . . . , Nr, and at destination 
D, respectively. We note that Q needs CSI about the source-to-relay and the relay-to-destination channels 
at relay and destination nodes, respectively. In this paper, we assume that CSI is perfectly known at the 
receiver while it is not known at the transmitter. This is obtained through adequate training [22]. 

After estimating bstR^ and 63^0, the destination D keeps the demodulated bit for further processing, as 
described in Section |llll while the relays initiate the relaying phase. More specifically, the generic relay, 
Rq, performs the following three operations: i) it applies binary NC on the set of demodulated bits bstR^ 
for t = 1,2, Ns; ii) it remodulates the network-coded bit by using BPSK modulation; and iii) it 
transmits the modulated bit to the destination D during time-slot T/Vg+g for g = 1, 2, . . . , Nji. Once again, 
we emphasize that all the demodulated bits are considered in this phase, even though they are wrongly 
detected, i.e., 65^ 7^ bStRq- ^s NC is concerned, we denote the network-coded bit at relay Rq 

by = {bs,R,,bs^R^, . . . ,bs^^R^ = gsiR^SiR, ® gs2R,bs2R, © ••• ® 9s,~,sR,bsi,^R,, where: i) 
fjl^ (•) denotes the encoding function at relay Rq; ii) © denotes exclusive OR (XOR) operations; and iii) 
gR, = [9SiR,,9S2R,,- ■ ■,9Sr,sR,]'^ is *e bin^y encoding vector at relay Rq [4], where gs,R^ G {0, 1} for 
t = 1,2,..., Ns. From this notation, it follows that only a sub-set of received bits are actually network- 
coded at relay Rq, i.e., only those bits for which gs^R^ = 1 for t = 1, 2, . . . ,Ns. Thus, our system setup 
is very general: no assumptions are made on g^^ for (7 = 1, 2, . . . ,Nji, and the encoding functions /r^ (•) 

'The framework proposed in this paper is applicable to other fading distributions. However, to keep the analytical development more concise and 
focused, we consider Rayleigh fading only. In Appendix |l] we provide some comments on how to extend the analysis to other fading distributions. 
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can be different at each relay. The goal of this paper is to understand how a given choice of these functions 
affect the end-to-end performance, as well as to provide guidelines for their design and optimization. 
Thus, the signal received at destination D in time-slot T^^^q after NC and modulation is ((7 = 1, 2, ... , Nr): 

VRgD = hu^DXRq + n_R,D (3) 

where xr^ = \/Em (l — 2&ij,). Let us note that the average transmit energy of each relay node is the same 
as the average transmit energy of each source node, i.e., Em- This uniform energy-allocation scheme stems 
from the assumption of no CSI at the transmitter. The impact of optimal energy allocation is postponed 
to future research [25]. Thus, the total average transmit energy for broadcasting and relaying phases is 
Et = Em {Ns + Nji), while the average transmit energy per network node is Ea = Et/{Ns + Nr) = Em- 

III. Receiver Design 

In this section, we develop a demodulator at the destination D which is robust to the error propagation 
problem caused by forwarding wrong detected bits from the relays. As explained in [1, pp. 18-20], the main 
goal is to improve the end-to-end performance by jointly performing demodulation and network decoding. 
To this end, we exploit the ML-optimum approach, which is composed of two main steps. 

a) Step 1: Upon reception of i/r^d in time-slot T/Vg+g, the destination D computes: 

f>R,D = argmin < Ij/flgD - V^B^/iiigD (l - 26^^ ) I | (4) 

where bR^D is the estimate of bR^. Two important comments are worth being made. 1) At the end of 
broadcasting and relaying phases, the destination D has Ns + Nr estimated bits, i.e., bstO for t = 
1,2,..., Ns from Q and bR^D for g = 1, 2, . . . , Nr from which can be seen as hard-decision estimates 
of all the bits transmitted in the network. These estimates are exploited in Step 2, as described below, to 
retrieve the information bits emitted by the sources and by taking into account NC operations performed at 
the relays. 2) Hard-decision demodulation is performed before network decoding, but, as we will better show 
in Step 2 below, the demodulator will take into account the reliability of these estimates when performing 
network decoding. Similar to [1, pp. 18-20], we will show that this is instrumental to achieve full-diversity. 

b) Step 2: In this step we take advantage of physical layer methods to develop network demodulation 
schemes that are robust to the error propagation problem [1], [8], and, thus, to the injection into the network, 
according to Q and Q, of wrong demodulated bits. This demodulator can be seen as a generalization of 
diversity-achieving demodulators for cooperative networks without NC [16]. The reader can notice that the 
proposed approach belong to the family of channel-aware detectors [54], [55]. To the best of the authors 
knowledge, the only notable paper which has recently extended these decoders to cooperative networks with 
NC is [25]. However, a single relay node and a fixed network code are considered in [25]. 
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Using the ML criterion, D demodulates the bits 65^ {t = 1,2, . . . , Ns) of the Ns sources as [52]: 



arg max 
6Si6{0,l},...,6s„ 6(0,1} 



{7'(Ssi,^'S2.---.^S„s|^Sii 



,bs. 



(a) 
OC 



6sie{o,i},...,f>s„ge{o,i} 



Ns 

Yl-P {bStD\bst 
t=l 



9=1 



Ylv(bR^D\bs„bS2 



,...,bf 



(5) 



(6) 

oc 



i,Sje{0,i},...,6sjv^ e{o,i} (^t=i 



arg max S Xl^'^ (^StD Ust)) + XI ( bfl,D Usi , ^Ss . ■ • ■ .^SNg)) 



where: i) (X| y) denotes the conditional Probability Density Function (PDF) of RV X given RV m ii) oc 

(a) 

denotes "proportional to"; iii) oc is obtained from the Bayes theorem, by exploiting the independence of the 

(b) 

detection events in each time-slot, and by taking into account that the emitted bits are equiprobable; and iv) oc 
is obtained by moving to the logarithm domain, which preserves optimality. Due to NC operations, in the sec- 
ond summation in the third row of ^ each addend is conditioned upon all the bits emitted from the source. In 



particular, from Section |II1 we have: V ( Sr^d 



bsi,bs2 



V ban 



fRg (bsi,bs2, ■ ■ ■ , bst. 



bR^j with bn^ = Jr^ {bsi,bs2,' 
The conditional probabilities in ([5]) can be computed as follows. By direct inspection, it follows that bstD 
for t = 1, 2, . . . , Ns turns out to be the outcome of a Binary Symmetric Channel (BSC) with cross-over prob- 
ability Ps,D = Pr {bs,D ^ bs,] = Q ^ 2{Em/ NQ)\hs,D?^ , where Q (x) = (l/^/2^) /+°" exp [-1^2) dt 
is the Q-function, Pr {•} denotes probability, and the last equality is due to using BPSK modulation. Accord- 



ingly, V[bs,D 



bst ) follows a Bernoulli distribution, i.e., V I bs^D 



55, 



{^-Ps^d) 



bstD-bst\ A^StD-bs 
^StD 



Similar arguments can be used to compute V [ bR d 



bR^ j . In particular, for q = 1,2,..., Nr, we have 



V bRn 



bR. 



l-Ps 1 : JVs D 



,1 



-bR, 



Si n R d ■ However, in this case the cross-over probabil- 



ity Psi:NsR,D = Pr |6/?,D / fR, {bs, , 652 , • • • , &5ivs ) } is no longer related to a single-hop link, but it must 
be computed by taking into account: i) dual-hop DemF protocol; and ii) NC operations performed at each 
relay node. To emphasize this fact, we use the subscript Si:N^RqD, where Si-n^ is a short-hand to denote 
the Ns sources of the network. This probability is better defined and computed in Section IIII-AI 

bR ] in (111), the ML-optimum demodulator simplifies, after 



By substituting V [bstD 



bs, and P bR^o 



some algebra and by neglecting some terms that have no effect on the demodulation metric, as: 

( Ns 

|6si>6s2.--->^Sjv<, « argmin ['^StD \bstD - bst\ 

bsie{0.i},...,is„ £{0,1} I t=i 



Nr 
9 = 1 



R„D lORnD 



(6) 



where ws,D = ln[{l - Ps,D)/PstD] and ws,,n^r^d = In [{l - PslnsR^d) / Ps^:MsR,D] for t = 1,2, Ns 
and q = 1,2, . . . , Nr, respectively. 

Three comments about Q are worth being made. 1) We can notice an evident resemblance with the 
well-known Chase combiner [48, Eq. (13)]. In spite of the similar structure, two fundamental differences 
exist between the original Chase combiner and i) the Chase combiner does not consider dual-hop 



^Throughout this paper, the PDF of RV X given RV Y is denoted either by V (X\ Y) or by Vx {■{ Y)- 
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networks, which means that all the packets reach the destination through direct links; and ii) the effect of 
error propagation caused by relaying and NC is not considered in the Chase combiner. These two differences 
are very important for two reasons: i) the detector in Q needs more CSI to work properly; and ii) the end- 
to-end performance of Q is affected by relaying and NC operations. Thus, the analysis of the performance 
of ^ requires new analytical methodologies, as we will better describe in Section |IVl 2) For large Ns and 
Nr, the complexity of Q can be quite involving. As suggested in [1, p. 19], this issue can be mitigated 
by using near-optimum demodulation methods {e.g., sphere decoding [56]), which attain ML optimality 
with an affordable complexity. 3) The demodulator in Q needs closed-form expressions of the cross-over 
probabilities Psi i^^r^d, which, in Section Ull-Ai is shown to depend on the CSI of the source-to-relay 
links, and on the NC operations performed at the relay nodes. In general, the estimation of this CSI requires 
some overhead [22]. In Section |Vl we will analyze the impact of CSI on the achievable diversity order. 

A. Cross— Over Probabilities of DemF— based Dual— Hop Networks with Binary NC 

In this section, Psi i^^r^d is computed in closed-form. Propo^/fton [7] summarizes the main result. 

Proposition 1: Let us consider system model and notation in Section |II] and Section |llll The exact cross- 
over probabiUty, Psi^N^iJ,!), for arbitrary binary encoding vectors is Ps^.j^^r.d = Psi..nsR<, + ^h,d - 
2Ps,.,,n^PR^D, where Pr^d = Q f a/2 (^^/iVo) |/ir,d|'V Ps,r, = Q ( \/2 (Em/ No) \hs,Rf] , and: 



(7) 



Ns 

t=l 



Ns 
r=t + l 



Proof: For the generic relay Rg, the end-to-end system can be seen as a dual-hop network where: i) 
the first hop is given by an equivalent wireless link with 6^^'' = Jr^ [bs^ ,bs2, ■ ■ ■ , &Sjvs ) its input and 
bRg = fRq (bsiR,,bs2R,,- ■ ■ , bsN^R,^ at its output, respectively; ii) the second hop is given by the wireless 
link with = /r^ (bsiR,,bs2R,, ■ ■ • at its input and Ir^d in © at its output. Thus, Ps^ ,^^r,d 

is given by Ps^,r,^R,D = Pr jftij^D / ^i^^"*}' which, by using [57, Eq. 23], is equal to: 

Ps,..^^R,D = Prjbfl, +Pr{S«,„ ^ b«J - 2Pr{bH, ^ fe^™)} Pr{6B,D 7^ ^fl, } (8) 

In ©, Fr^bR^D^bR^} = Pr^d = Q (^2{Em/NQ) as it is the error probabiUty of a 

single-hop link. On the other hand, Pr ^r^ 7^ ~ Psi NgRq explicitly written as: 

PSi,T,gRq = PrjaSiflg^Siii, ©952^5 ^Sa-B, © • • ■ © SSjvs-Bg^Sjvgfl, QSiR^bs^R^ © c/Ss-R^^Sa-R, © • • ■ © SSjvs -R, ''S„s B, } (9) 



Let us now introduce the notation (t = 1, 2, . . . , Ns)'- 

(PSi.,tRq = ^'^{sSirJ'SiR, ®9S2RgbS2Rg ® 9StRgbstRg + 9SiRgbsiRg ® 9S2Rgbs2Rg © ■ • ■ © 9St if., f'St } 

^StRg'^''^ = ^^{aStRqbStRq 9StRqbstRg} = 9StRq Prj^Stfl, ¥= bstRg} = 9StRqPStRg 

With = P^r"'^ = if * = 1- Furthermore, similar to Pr^d, Ps^r, = Pr {bs,R, + ft^.i?,} 
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PstRg = Q 2 {Em/ Nq) \ hs^R^ I J . By taking into account the properties of the XOR operator, (O can 
be computed by using the following chain of recurrence relations: 




A closed-form solution of a recurrence relation similar to ([TT]) has recently been given in [58] for multi- 
hop networks. In particular, by using [58, Eq. (9)], (|7]) can be obtained. This concludes the proof. □ 

Proposition [7] is instrumental for an efficient implementation of Q. Furthermore, the proof sheds lights 
on the fundamental behavior of NC over fading channels. In fact, by comparing ^ and [58, Eq. (9)], we 
notice that the cumulative error due to performing NC on wrong demodulated bits at the relay is equivalent 
to the error propagation problem in multi-hop networks. In other words, if the relay performs NC on the data 
received from 1 < A'^^ < Ns sources, then the error probability of the network-coded data is the same as a 
multi-hop network with A'^^ hops having fading channels given by the source-to-relay links. When adding 
the relay-to-destination link, the end-to-end network behaves like a Ng + 1 multi-hop network. In other 
words, Propo5/rio« [7] clearly states that the larger the number of network-coded sources is {i.e., the larger 
the number of non-zero elements of the encoding vector g/j^), the more important the error propagation 
effect might be. In summary, Pro;?o5jY/o?i [7] provides a simple, compact, and intuitive characterization of the 
error propagation caused by DemF relaying and NC over fading channels. 

IV. Performance Analysis - ABEP 

In this section, we provide closed-form expressions of the ABEP for each source of the network. The 
framework takes into account the DemF relay protocol and the characteristics of the network code. The 
departing point of our analysis consists in recognizing that, according to Section |lll the network code can 
be seen as a {Ns + NR)-\ong distributed linear block code, whose first Ns bits can be seen as systematic 
information bits, and the last Nr bits can be seen as the parity (redundant) bits. However, there are two 
fundamental differences between the system model under analysis and classical linear block codes [52], 
[59]: i) the system model in Section HIl encompasses a dual-hop network, while state-of-the-art analysis 
of classical codes usually considers single-hop transmission; and ii) coding is not performed at the source 
nodes, but it is performed at the relay nodes. Due to the distributed nature of the network code and the 
assumption of realistic fading channels, encoding operations at the relays are inherently error-prone, as shown 
in Proposition [7] Accordingly, new frameworks are needed to characterize the end-to-end performance of 
dual-hop networks with NC, similar to the many frameworks that have been developed for cooperative/multi- 
hop networks without NC [58], [60]. Note that the frameworks proposed in [44] and [45] are not applicable 
to our setup since fading is neglected and no diversity-achieving demodulators are investigated. 
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Owing to the inherent similarly between the system model in Section |ll] and distributed linear block codes, 
we use union-bound methods to compute the ABEP [53, Eq. (12.44)]. The main difference with respect 
to state-of-the-art frameworks is the computation of each individual Average Pairwise Error Probability 
(APEP), which must account for the DemF protocol and for the error propagation introduced by NC. 
Furthermore, in this paper we are interested in computing the ABEP of each source of the network instead 
of considering frame or codeword error probabilities, as it is usually done for linear block codes [52], [53]. 
The reason is that in our distributed system each source transmits independent information flows, and we 
are interested in characterizing the error performance of each of them. 

Using the union-bound for equiprobable transmitted bits [53, Eq. (12.44)], the ABEP of source St is: 

ABEPs,<^X: E E E E E [APEP(c^c)A(c[t],c[t])] 

- i E E [APEP (c ^ c) A (c [t] , c m 

b b 

where: i) = is a short-hand to avoid multi-fold summations; ii) (•) denotes transpose operations; iii) Inxn 
is an n X n identity matrix; iv) b = [65^,65,, . . .,bs„^]'^ and b = [bs,,bs^,. ■ ■ ,bs„^]'^; v) c = G0b and 
c = G b, where "0" indicates that matrix operations (additions and multiplications) are performed in the 
Galois field GF(2), G = [gRi , SR2 ) ■ ■ ■ ; SRn ] ^ t^e Nr x A'^s' matrix containing the encoding vectors of 
all the relays, and G = lArgxAfslG is the (Ns + Nr) x Ns generator matrix of the whole distributed 
network code; vi) z [m] is the m-th entry of vector z; vii) A = 1 — A where A (•, •) is the 

Kronecker delta function, i.e., A (x, y) = 1 if x = y and A (x, y) = elsewhere; and viii) APEP (c — c) is 
the probability, averaged over fading channel statistics, of detecting c when, instead, c is actually transmitted, 
and these are the only two codewords possibly being transmitted. The Kronecker delta function takes into 
account that a wrong demodulated codeword might not result in an error for the source, St, under analysis. 

The next step is the computation of the APEP for a generic pair of distributed codewords. We proceed 
in two steps: i) the PEP conditioned on fading channels is computed; and ii) the conditioning is removed. 

A. Computation of PEP (c — )■ c) 

The decision metric in Q can be rewritten in a more compact form as follows: 

E ('"StD pSti3 - ^St|) + E ("'Si:Ns-R<,-D pflgD - fcfl, I) = ^ {w [m] |c[m] - c [m]|} = A(c) (13) 
t—l Q=l m=l 

where we have defined: w = [ws^d^ws^d, ■ ■ ■ ,ws,,^d,ws^.,,,^r^d,ws^,^^R2D, ■ ■ ■ ,ws^,^^r^^dY , b 



T _ 



T 



bSiD,bs2D, • • • , bsj^^D,bR^D, bR^D, ■■■ , bR^^D , b = bs^D,bs^Di ■■■ , bsj^^D,bR^D, bR^D, ■■■ , bR^^D 

c = Gb, and c = Gb. From the PEP i.e., PEP (c ^ c) = Pr {A (c) > A (c)}, is: 

(Ns+Nn ] la) i 1 

PEP (c c) = Pr ^ w[m] (|cH - cHI - |c[m] - c[m]|) > > = Pr ^ w [m] (|c [m] - c [m] | - |c [m] - c [m] |) > 

[ ™=i J [mee{c,c) J 

(14) 
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where == is obtained by taking into account that D (c [m] , c [m]) = w [m] (|c [m] — c [m] \ — |c [m] — c [m]\) 
contributes to the summation if and only if c [m] 7^ c[m], and, thus, the summation considers only the 
elements in the set 0(c, c) = {m|c[m] 7^c[m]}. The cardinality, card{-}, of 0(c,c) is given by the 



Ns+Nr 
m=l 



c \m\ — c \m\ 



Hamming distance between c and c, i.e., card {6 (c, c)} = dn (c, c) = ^ 

By conditioning on m G (c, c), it can be shown, by direct inspection, that D (c [m] , c [m]) is a discrete 
RV which can only assume values w [m] and —w [m] with probability P [m] and 1 — P [m], respectively, 
where P = [Ps,d, PS2D, ■ ■ ■ , Ps^sD, Ps^.m^RiD, Ps,.nsR2D, ■ ■ ■ , Psi..m^Rm^d] ^ is the vector of cross-over 
probabilities computed in Section Hill Accordingly, the conditional PDF of RV D (c [m] , c [m]) is: 



7'D(c[™],c[m]) (?| m e e (c, c)) = P [m] 5 - w [m]) + (1 - P [m]) 5 + w [rn]) 



(15) 



where 6 (•) denotes the Dirac delta function. Since the RVs D (c [m] , c [m]) are independent for m G © (c, c), 
then D (c, c) = X]mee(c c) ^ t^l ' ^ t"^]) ^ given by the convolution of the dn (c, c) PDFs of the 
individual RVs D (c [m] , c [m]). More specifically, let us denote by ni(c,c) = |"^(c^c)' "^(c^c)' • • • ' "^(c e)^'*^^^ } 
the specific set of dn (c, c) indexes such that m G © (c, c). Then, the PDF of D (c, c) can be written as: 



V 



D(c,c) 



(0 



(16) 



where (8) denotes convolution operations. Thus, by definition, the PEP in ([T4l ) can be computed as PEP (c — )• c) = 
j+co -p^^^ (^^^ closed-form expression of this PEP is given in Proposition |2l 

Proposition 2: Let Pd(c,c) (•) in Gil, for high-SNR {i.e., Em/No oo), PEP (c c) = Jq°° Vd{c,c) (0 
can be tightly upper-bounded as follows: 



J„(o,c) 

n ^ 




fcG*i(c,c 


k = l 


mill ■ 




E 

fc=i 

):e<l>3(c,c) 


mill • 





h=l 
h^k 



dff (c,c) 

n 

h = l 

ft^Vfc[l],)i^Vfc[2] 



,(-fc[21)l 
(c.c) 



(t:,c) 



te4.3(o,c) 



djj(c,c) 

n 

h=\ 



(17) 



lLdH{=,e)/2j^ 

'=«^*L<1h(<=Tc)/2J<<='^) 



(c,c) 



^(-fc[LrfH(<=.e)/2j])l 

(c,c) 



n phs:-'.,] 

'•5^^'fc[^'l,'lVvfc[2] 

[dj,(c,c)/2] 



where: i) ( ) is the binomial coefficient; ii) [-J is the floor integer part; iii) '^n (c, c) is a set of in- 
dexes defined as (c, c) = \^k\k + Y^lzX{'^"^^'^'^) < 2"'«(^'^)-i - l|; and iv) ^>„(c,c) is the set of 
all possible combinations of the indexes in ni(c,c) taken in sets of n, and it is defined as (c, c) = 
{vfcl Vfc = {vfc [1] , Vjt [2] , . . . , Vjt [n]}}, where is its /c-th element, i.e., the fc-th combination of the 
indexes in m((.,c)- The cardinality of <!>„ (c,c) is card{<I>„ (c, c)} = (^"^'^'*)- 
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Proof: We proceed in two steps: i) first, we describe the step-by-step methodology to compute PEP (c — )• c) 
in (fTTl) for du (c, c) = 3; and ii) then, we describe how the approach can be generalized to generic du (c, c). 



Let us start with du (c, c) = 3. In this case, we have m 



(c,c) 



,(1) 



,(2) 



,(3) 



}■ 



and V 



Dice 



(•) in 



([T6l ) can be computed by using some properties of the Dirac delta function. By doing so, and substituting 
the obtained PDF in PEP (c ^ c) = f^°^ 'Poicc) (0 we get: 



PEP (c -> c) = P 






















,])( 


L - P 


Nc:.,])(- 






+ P 


HZ 


,](^ 








s.,])«(wHc:., 


-Ky-wH^,^,]) 


+ P 


I - (2) 
(o,c 


)](^ 




HZ 








+ P 


r - (3) 

L (c,c 


)](^ 






])(.-p[. 




,]-wHcy+wH2y) 


+ P 
















+ P 








.](^ 








+ P 


r - (2) 

L (o,t 






.](^ 


--He!.,]) 





(18) 



where 71 (x) = Jq °° ^ — x) d^ is the Heaviside function: ?^ (x) = 1 if x > and ^ (x) = elsewhere. 

The PEP in ([TSl l can be simplified and can be written in a form that is more useful to compute the average 
over fading statistics. The main considerations to this end are as follows: i) since, by definition (see (|6]l). 



w [m] > for m = 1, 2, . . . , d// (c, c), then Ti { — ^ 



iff (c,c) 



W 



_ (fc) 
(c,c) 



Oand?^ 



^dff (c,c) 



W 



m 



(k) 

{c,c) 



1 for any dn (c, c) and for any m(c g); and ii) in the high-SNR regime, the BEP in (1181 ) can be tightly 
upper-bounded by recognizing that 1 — P [m] — )• 1 for m = 1, 2, . . . , d// (c, c). Furthermore, by exploiting 
i) and ii), the resulting terms containing the Heaviside function can be grouped in three pairs of two addends 
each. For example, a pair in ([TSl l is Z = Zi + Z2 with: 



\ Zi=P 


{c,c) 


[ Z2=F 


■ _ (3) 
{c,c) 



(1) 

c,c) 



v. 



{c,c) 



- 3 
(c,c) 



(19) 



(c,c)J I (c,c)J I (c,c)J 

while the other two pairs can be obtained by direct inspection of ([TSll accordingly. 
For generic du (c, c), pairs as shown in (l20b can be obtained: 





= n ■ 


P 




l«( 












Z2 


= n • 


P 




|«( 




keA 









fee.4 



■ (k) 
{c,c) 



w \m 



k£A 



■ (fc) 

{c,c) 



(20) 



where A and ^ are two sets of indexes such that m^f, = ^ U ^ and ^ n ^ = 0. 

By taking into account that, for high-SNR, we have w [m] = In [(1 — P [m])/P 
and from the definition of Heaviside function, ^ (•), Zi and Z2 in (l20l ) simplify as: 



m\] 



■In (P [m]). 



Zi 



n pHc%) 



if n p [™['!.)] > n p ['7^1'!.) 

keA L ^ ' j.g_4 L \ > / 
elsewhere 



Z2 



n p [ml^Ul if n p 







elsewhere 



(fe) 

(c,c) 



< n p 

keA 



{c,E)J 



(21) 
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Thus, for high-SNR, Z = Zi + Z2 can be re-written as follows: 

n if n p[<.y < n_p[ 
if n p [<y < n p 



n P |m 

keA 



(fc) 

(c,c) 



n ] 

fce^ 



mm < P 

keA 



(c,c)J ' i 1_ [ (c.c) 
keA 



(22) 



In conclusion, by exploiting the properties of the Heaviside function, (•) and the high-SNR approxi- 
mation in (I2T] ). the PEP in (ITSl ) can be tightly upper-bounded as follows: 



PEP (c ^ c) ^ P m 



(c,c) 



P m 



- (2) 
{c,c) 



P m 



- (3) 
\c,c) 



j + min I 



P \rh 



(CC) 



,P m 



=,(2) 



(cc) 

P |m 



(1) 

(c,c) 



(23) 



The result in (1231 ) represents the first part of our proof, and allows us to explain two main aspects of (ITtI ): 
i) its validity and accuracy for high-SNRs only, as some approximations are used; and ii) the presence of the 
min{-, •} function, which comes from grouping pairs of addends, and by exploiting definition and properties 
of the Heaviside function. The second step is to provide a justification of ([TT] ) for arbitrary du (c, c). First, let 
us emphasize that, when possible, the proof for (c, c) = 3 has been given for arbitrary dn (c, c), which 
provides a first sound proof of the generality of our approach. Second, we emphasize that the interested reader 
might repeat the same steps as for the case study with du (c, c) = 3 for arbitrary du (c, c) and eventually 
lead to ([TT] ). The only difficulty if the large number of terms arising when computing the convolution in 
([T6l ). So, here we provide only some guidelines to understand ([Tt] ). The first thing to observe is that (|23]) 
can be obtained from ([TT] ). and, more specifically, it is given by the first two addends in the right-hand side 
of ([TT]) . The other terms come from the fact that, for dn (c, c) > 3, in (l22l ) we have to consider all possible 



combinations of the indexes ni(c,c) taken in sets of 1, 2, 3, etc., since A and A in (1201 ) are a partition of the 
dn (c, c) indexes in m(c,c)- This explains the presence of all the other summations in ([TT] ). along with the 
upper limit of each of them. The reason why the upper limit of the last summation is (|^£;^^c'eV2j) '^^^ 



to the equality min|]^^g^P 



(c,c) 



keA- 



m 



(k) 

(c,c) 



} = min{n,g^P He,y 'H^e^P Hey}' 



and because only one of these latter terms is explicitly present in ([T6l ). Furthermore, the need to compute 
all possible combinations of the indexes m(c,c) clearly explains the definition of (c, c) in ([TT] ). The 
only thing left is to understand why in each summation the index k must belong to the set ^„ (c, c). The 
motivation is as follows. When computing the convolution in (fT6l ). the total number of addends in the final 
result is 2'^"(^'^). In fact, the convolution of 2'^"^'^''^^ PDFs is computed, each one given by the summation of 



two terms. Among all these 2'^"^'^''^) terms, 11^=1 
separately in (fTTl ). More specifically, fl 

while ntf f 1 - P 



dff(c,c) 



dff (c,c) ■ 
k=l 



- (fe) 
(c,c) 



and nt"JrUi-F 



m 



(k) 

(c,c) 



are treated 



(k) 

(cc) 



(c,c)^ 



is explicitly shown in ([TT] ) as the first addend. 



in zero because of the properties of the Heaviside function. The remaining 
2rfff(c,c) _2 grouped in pairs of two addends, as shown in (l20l) . Furthermore, each pair reduces to only one 
addend as shown in (l22l ). Accordingly, the number of terms in (fTTl ) cannot be larger than {2'^h{c,c) _ 2)/2 = 
2rfff(c,c)-i _ Y other words, when the cumulative inequality in (c, c) is no longer satisfied, we can 
stop computing the summations in (fTTl) . This concludes the proof. □ 
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Proposition\2\is very general and can be applied to any (c, c). However, it is not an exact result, as it 
holds for high-SNR only. For the special case dn (c, c) = 2, an exact expression of the PEP in ([141 ) can be 
obtained, which, in general, has to be preferred as it is accurate for any SNR. In Corollary |7] we provide 
the exact expression of the PEP in ([141 ) without any high-SNR approximations. 

Corollary 1: If dn (c, c) = 2, then m(c,c) = |"^(c^c)' "^(c^c)} ^^'^ (1141 ) is equal to: 

PEP (c c) = min {P [rhfj-^] , P [mg_j] } (24) 

Proof: The proof follows from analytical steps similar to ( fTSl ) in Proposition \2\ In particular, we have: 



PEP (c ^ c) = P P [rn[%^] + P [rn^]}.^] (l - P [rn'^-^^])n (w [fh^^].^^] - w 



(25) 



UnUke Proposition \2\ there is no need to exploit the high-SNR approximation 1 — P [m] — 1. On the 
contrary, by using the properties of the Heaviside function, ((2T1) . and (|22l) . we get: 

PEP (c ^ .) - I ^ ['^('^)] ^ i'^^-)] ^ - ^ ['^('^)]) ^ i'^'-)] ^ ^ ['^('^)] ^ ['^(-)] < ^ ['^C'^)] (26) 
|p P [^|^y+(l-P P[<y=p[<^^^ ifp[^«,]<p[^|^y 

which clearly leads to (|24l) . This concludes the proof. □ 
We note that the main difference between (ITtI ) and (l24l ) is the absence in (l24b of the first addend in (fTTl ). 
In fact, this addend simplifies if the high-SNR approximation 1 — P [m] — )• 1 is not used in (l26t . This 
provides a better (and exact) estimate of the PEP. However, this procedure cannot be readily generalized to 
network codes with dn (c, c) > 3, without having a more complicated expression of the PEP, which is not 
useful for further analysis, and, more specifically, to remove the conditioning over fading statistics. 

B. Computation of APEP (c — > c) 

The aim of this section is to provide a closed-form and insightful expression of the APEP, i.e., to average 
the PEP in (fTTl ) over fading channel statistics. In spite of the apparent complexity of (ITtI) . Proposition \3\ 
shows that a surprisingly simple, compact, and insightful result can be obtained for i.n.i.d. fading. 

Proposition 3: Let us consider the Rayleigh fading channel model introduced in Section Hi] The APEP, 
APEP (c ^ c) = Eh {PEP (c c)}, is as follows: 



APEP (c ^ c) ^ 4 

^ ' ^ NoJ 



K \ -dff (c,c) 



1 + 20fr Uh (c, c) + - J2 



^ Ldff(c,E)/2J j^{dH(c,c)) 



d 



l+\)V{dH (c,c)-d+i) 



Ns + Nr 

X 

m— 1 



(27) 



1 T • i^ 



where: 



d 



2dff (c,c)-l - 1 - J] ^dij{c,c)-j J2 (dH(c,c)^ ^ 2'*ff(<='<=)-l - 1 

e — 1 e — 1 

and: i) Eh{-} denotes the expectation operator computed over all fading gains of the network model 
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introduced in Section HIl ii) x {C} = 1 if = and x {C} = ? if C / 0; iii) Ac,c = c c = (G b) 



(G0b); iv) S^^i^ = Ssd + Srd + s(^^; v) Ssd = [l/^z?, • - 1 A 



T 



[Sgj-,, OixAT^j] , where Oixn is a 1 x n all-zero vector; vi) Ssr^ 



-'SR 



T 



.(G) 
IxNs, [ ^SR 



and Ssd 



; and vii) Srd 



, and S^p^ = 

xA^s'^rd] • Finally, we emphasize that in S^^^ 
usual matrix operations are used and arithmetic is not in GF(2). 

Proof: From the definition of APEP, i.e., APEP (c — ^ c) = Eh {PEP (c — ^ c)} and the linearity property 
of the expectation operator, it follows that two types of terms in dTVl ) have to be analyzed: 



rfff {c,c) 



h < mm 



■ {*:) 
{c,c) 



(c,c) 



(29) 



An asymptotically-tight (for Em/No — )• oo) approximation of Ti and T2 in ( |29l ) can be obtained by using 
Lemma [7] and Lemma\2\m Appendix H In particular, for high-SNR, ( |29l ) simplifies as follows: 



dH-(c,c) Ns+Na 




2v^r(d„(c,c)+i) 

r(d+i)r(d„{c,c)-d+i) 



n X 

m— 1 



{Ac 



[ml S 



(<3) 



m} 



(30) 



where d = card {A} denotes the cardinality of set A. 

From (l30l ). equation (|27] ) can be obtained from ([TT] ) as follows: i) the first addend in (ITtI) is Ti in (l29l) 
and, thus, it can directly be obtained from (l30l ): ii) each min {•, •} term in (fTTl ) corresponds to r2 in ( [29l ) and, 
thus, it can directly be obtained from (l30l) : and iii) by carefully studying T2 in (l30l) . it can be noticed that it is 
independent of the particular sub-set of indexes in A and A, as defined in ( |29l ). The only thing which matters 
is the number of indexes in A and in A, i.e., their cardinality card {A} = d and card {A} = dn (c, c) — d. 



respectively. For example, if du (c, c) = 3 in ( |23] ). then Eh |min 

}} = Eh {min {P 



Eh { min < P 



(2) 
(c,c) 



(1) 
(c,c) 



(3) 
(c,c) 



m 



- (1) 

(3) 
(c,c) 



- (2) 
(c,c) 

(1) 

(c,c) 



m 



m 



(2) 
(Cc) 



(3) 
(c,c) 



}}■ 



}} 



This 



remark holds for generic i.n.i.d. channels, and it implies the identity (for n = 1, 2, . . . , [dn (c, c) /2j): 



Eh < 



E 



1 

fcG*7i(c,c) 

vfce«„(c,c) 



mm < f3 P 



{c,c) 



n p ['^fe'; 



(31) 



where r2 is given in (Unil, and AT 



(dff (c,c)) 



is the number of terms in ( [281 ) that are actually summed in ( [3T| ). 



By putting together these considerations, and by taking into account that there are [dn (c, c) /2j summa- 
tions with different card {A} = din (fTTl l. we obtain (|27] ). The only missing thing in our proof is to show that 
j^idnicc)) closed-form expression given in ( [28] ). This result follows from the definition of (c, c) 

for n = 1, 2, . . . , [dn (c, c) /2j in In fact, since ^„ (c, c) = | A:| A: + Y^lZl {'^"t^^) < 2'^"^"^^^-^ - l}, 
the number of elements in each summation in (|3TI ) is: i) either {^"^^''^^) , if we have not reached the maximum 
number of indexes that can be summed, i.e., 2'^«(*^'*^)~^ — 1; ii) or, in the last summation, the remaining 
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indexes if the cumulative summation in ^„ (c, c) exceeds this maximum number of indexes. Equation (1281) 
summarizes in formulas these two cases. This concludes the proof. □ 
Similar to Proposition |2j the exact APEP can be obtained if cLh (c, c) = 2, as given in Corollary |2l 
Corollary 2: Let us consider the Rayleigh fading channel model introduced in Section |lll Then, APEP (c — ) 
Eh {PEP (c c)} with PEP (c c) given in for du (c, c) = 2 is as follows: 

APEP(c-^c)= (V2^j n x{Ac,eMs(^1,m} (32) 

^ m— 1 

where the same symbols and notation as in Proposition \3\ are used. 

Proof: It follows from (|27] ) with dn (c, c) = 2, by neglecting the "1" term as shown in Corollary \l] □ 

C. Particular Fading Channels 

Proposition \3\ is general and it can be applied to arbitrary i.n.i.d fading channels and network topologies 
with generic binary NC. However, it is interesting to see what happens to the network performance for some 
special channel models and operating conditions, which are often studied to shed lights on the fundamental 
behavior of complex systems. In this section, we are interested in providing some simplified results for 
three notable scenarios of interest: i) i.i.d. fading, where we have = (Tq for every wireless link; ii) 
i.n.i.d. fading with high-reliable source-to-relay Unks, which is often assumed to simplify the analysis, but, 
as described in Proposition I///-AI it does not account for the eiTor propagation effect due to NC; and iii) 
i.i.d. scenario with high-reliable source-to-relay links. The end-to-end APEP of these three scenarios is 
summarized in Corollary \3\ Corollary and Corollary \5\ respectively. 

Corollary 3: If the fading channels are i.i.d. with a^y — <^0' '^^^'^ APEP in Proposition \3\ and in 
Corollary^ can be simplified by taking into account the following identity: 

Ns+Nr _ Ns+Nr 

n x{Ac,.Ms(°)3H} = (-g)"'"'^''' n x{Ac,.Mg(o'M} (33) 

7n=l m—l 

where: i) U^Ns is a 1 x iV^ all-one vector; ii) g^^ = 1 + ^^J^ 9s,R, = 1 + N^'^^^ for q = 1,2, Nr, 
where N^^'^"'' is the number of sources whose data is network-coded at relay node Rg-, and iii) g^^^ = 

"i (0) (0) (0) 1^ 

Proof: It follows from Proposition \3\ with a^y — '^o- particular, 5]sD and 5]rd simplify to all-one 
vectors multiplied by 1/(Tq, and each entry of 5]g^^ reduces to the summation of the elements of the binary 
encoding vector used at each relay, which is equal to the number of network-coded sources. □ 

The result in ( [33l ) is very interesting as it clearly shows, through ]^'^^'^'^'*\ that the larger the number of 
network-coded sources is, the more pronounced the error propagation problem might be. Thus, depending 
on the quality of the fading channels, it might be more or less convenient to mix at each relay the data 
packets transmitted from all the sources. Further comments are postponed to Section IV-AI 

Corollary 4: Let us assume that the source-to-relay channels are very reliable, i.e., no demodulation errors 
at the relays. For example, this can be achieved either by using very powerful error correction codes on the 
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source-to-relay links, or when the relays are located very close to the sources. Then, 5]gj^J-, in Proposition\3\ 
and Corollary |2] simplifies to E^^j^ = Ssd + IIrd = l/ cts.d^ ■■■A/ ^^s^.d^ l/ '^r.d^ ■■■A/ '^r^^d ■ 

Proof: If the source-to-relay channels are very reliable, we have cr|^^ — )• oo for t = 1, 2, . . . , Ns and 
g = 1, 2, . . . , A'^/j. Thus, by definition, 5]gj^ — )• 0((7Vs+ArH)xi)- So, the simplified expression of 5]gj^^ follows 
by taking into account the definition of Ssd and Srd as block matrices. This concludes the proof. □ 

Two important conclusions can be drawn from Corollary HI First, we notice that the APEP is affected 
by the encoding operations performed at the relays only through the codeword's distance dn (c, c), which 
is the number of distinct elements between c and c. This provides a very simple criterion to choose the 
network code for performance optimization. Second, since 5]gp{^j-) [m] < 5]gj^j-, [m] for 

171 = 1,2,..., {Ns + Nr), then APEP (c c)| , < APEP (c 'c)| 2 which is an expected 

result, and it confirms that, to limit the error propagation due to NC operations, the source-to-relay links 
should be as reliable as possible. Further comments are postponed to Section IV-AI 

Corollary 5: If the fading channels on the source-to-destination and relay-to-destination links are i.i.d. 
with cr^y = cTg, and the source-to-relay channels are very reliable with no decoding errors at the relays, 
then Proposition\3\md Corollary \2\ can be simplified by taking into account the identity: 



Ns+Nn 



(34) 



Proof: It follows from Corollary |?1 which for i.i.d. source-to-destination and relay-to-destination links 
jives Sg^jj = (l/o"o) l((Ars+Ar^)xi)- Since there are dn (c, c) non-zero terms in Ac,c, we get (l34l) . □ 



V. Analysis of Diversity Order and Coding Gain 

To better understand the performance of the cooperative network under analysis, and to clearly showcase 
the impact of the distributed network code on the end-to-end performance, in this section we study diversity 
order and coding gain according to the definition given in [61]. In particular, we are interested in re-writing 



the end-to-end ABEP in ^ as ABEP^^ 



-G 



(St) 



, where ci^*^ and G^f*^ are coding 



gain and diversity order of S'f for t = 1, 2, . . . , Ns, respectively. This result is summarized in Proposition 

m 

Proposition 4: Given the ABEP in ([T2l) and the APEP in ( [27] ). diversity order and coding gain of St are: 



G 



(St) 



SV [t] 



1 

2™S 



E 

b.b 

dff(c,c) = SV[t] 



Ld„(c,c)/2J 



1 + 20Fr {da (c, c) + i) E 



r(d+i)r(d„(c,c)-d+i) 



Ns+Na 



n 



SRD I 



,]}) A(cW,cM) 



(35) 



where SV is known, in coding theory, as "Separation Vector" (SV) [63, Def. 1], and, for a given codebook 
C = {c| c = Gb, Vb}, its t-th entry, i.e., SV [t], is defined as the minimum Hamming distance between 
any pair of codewords c = Gb G C and c = Gb G C with different t-th bit, i.e., with b [t] / b [t]. 



18 



Proof: First of all, let us study G^f*\ From (l27l ) in Proposition\3\'we notice that the APEP has diversity 
order dn (c,c) [61], which is the Hamming distance between the pair of codewords c and c. Furthermore, 
from ([T2I ) we know that this APEP contributes to the ABEP of source St if and only if the t-th bits of c 
and c are different, i.e., if and only if c[t] 7^ c [t]. Since the network codes studied in this paper can be 
seen as systematic linear block codes, as explained in Section |II1 the latter condition implies b [t] 7^ b [t]. 
Accordingly, in ( |27l ) only the APEPs having a diversity order, i.e., a Hamming distance, equal to: 

d^J^'"' (t) = I (c, c)| cIh (c, c) < dH (c' , c') Vc, c, c' , c' eC with c[t]=/=c [t] and c' [t] c' [t] } (36) 

will dominate the performance for high-SNR. In fact, all the other APEPs will decay much faster with the 
SNR, thus providing a negligible contribution. In formulas, the ABEP in (fT2l ) can be re-written as: 

ABEPst < ^^X^[APEP(c->c)A(cM,c[t])] ^ ^ J2 [APEP (c c) A (c [t] , c [t])] (37) 

b b b,b 

d„(c,c) = d^'">(t) 

From dMJl and by definition [63, Def. 1], d^^''''^ {t) is exactly the i-th entry of SV, i.e., d^'"^ (t) = 
SV [t]. Thus, we have proved that the end-to-end diversity order of source St is equal to its SV. This result 
showcases that, depending on the used network code, different sources in the network have, in general, 
different diversity orders. This observation has important applications, as described in Section IV-AI Finally, 
the coding gain, Gi^*\ can be obtained through algebraic manipulations by substituting (|27] ) in (|37] ). and 



r (s ) 

equating the resulting expression to ABEP5^ — {Em/No) Gc 



This concludes the proof. □ 



A. Insights from the Analytical Framework 

Even though the overall analytical derivation and proof to get ( |271 ) in Proposition \3\ are quite analytically 
involving, the final expression of the APEP turns out to be very compact, elegant, and simple to compute. 
In Section [VTl via Monte Carlo simulations, we will substantiate its accuracy for high-SNR. In addition, 
the framework is very insightful, as it provides, via direct inspection, important considerations on how the 
network code affects the performance of the cooperative network, as well as how it can be optimized to 
improve the end-to-end performance. Important insights from the analytical framework are as follows. 

• End-to-end diversity order. As far as diversity is concerned, in Proposition |?] we have proved that each 
source can achieve a diversity order that is equal to the separation vector of the network code. This is a 
very important result as it shows that even though a dual-hop network is considered, which is prone to error 
propagation due to relaying and to demodulation errors that might happen at each relay node, the distance 
properties of the network code are still preserved as far as the end-to-end performance is concerned. This 
result allows us to conclude that, if we want to guarantee a given diversity order for a given source, we can 
use conventional linear block codes as network codes, and be sure that the end-to-end diversity order (and, 
thus, the error correction capabilities [49], [50], [63]) of these codes is preserved even in the presence of 
error propagation due to relaying and NC operations. This result and its proof is, to the best of the authors 
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knowledge, new, as it is often assumed a priori that the presence of error propagation does not affect the 
diversity properties of the network code [34]. On a practical point view, this result suggests that, as far as 
only the diversity order is concerned, the network codes can be designed by using the same optimization 
criteria as for single-hop networks. Finally, we note that the result obtained in this paper is more general 
than [18], as our proof is not based on the Singleton bound, and, more important, no ad hoc interleavers 
are needed to achieve a distributed diversity equal to the SV. 

• Comparison with single-hop network and classical coding theory. It is interesting to compare the result 
about the achievable diversity order in Proposition |?] with the diversity order that is achievable in single- 
hop networks. From [52, Sec. 14-6-1)], [53, Ch. 12], and [59, Sec. II], we know that single-hop networks 
operating in fully-interleaved fading channels and using soft-decision decoding have a diversity order that is 
equal to the minimum distance of the linear code. The result in Proposition^can be seen as a generalization 
of the analysis of single-hop networks in [52], [53], [59] to dual-hop networks with NC. It is important to 
emphasize that in our analysis we have taken into account realistic communication and channel conditions, 
which include demodulation errors at the relays and practical forwarding mechanisms. Also, our results are 
in agreement with [64], where the error correction properties of network codes for the single-source scenario 
have been studied, and a strong connection with classical coding theory has been established. Our analysis 
extends the analysis to multi-source networks, provides closed-form expressions of important performance 
metrics, and accounts for practical communication constraints. Finally, we note that even though relays and 
destination compute hard-decision estimates of the incoming signals and send them to the network-layer to 
exploit the redundancy introduced by cooperation and NC, the diversity order is the same as in single-hop 
networks with soft-decision decoding. The reason is that at the network-layer we take into account the 
reliability of each bit through a demodulator that resembles the Chase combiner [48] (see also Section [V-BI ). 

• Comparison with adaptive NC solutions. In Section IH we have mentioned that another class of network 
code designs aims at guaranteeing a given end-to-end diversity order without injecting erroneous packets 
into the network. In these solutions, the network code changes according to the detection outcome at each 
relay node. Results and analysis in [23], [27], and [28] have established a strong connection between the 
design of diversity-achieving network codes and linear block codes for erasure channels. More specifically, 
[23], [27], and [28] have shown that MDS codes can be used as network codes to achieve distributed diversity 
for erasure channels. The analysis conducted in the present paper complements design and optimization of 
network codes for erasure channels to the performance analysis and design of such codes for error channels, 
where all the bits are forwarded to the destination regardless of their reliability. 

• End-to-end coding gain. As far as the coding gain in Proposition^is concerned, and unlike the analysis 
of the diversity order, there are differences between single- and dual-hop networks with and without NC. 
In fact, in Corollary we have shown that both demodulation errors at the relays and dual-hop relaying 
introduce a coding gain loss if compared to single-hop transmissions. Thus, even though NC and relaying, via 
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a proper receiver design, do not reduce the diversity order inherently provided by the distributed network 
code, they do reduce the coding gain, which results in a performance degradation that depends on the 
quality of the source-to-relay channels. However, this performance degradation might be reduced, and even 
completely compensated, through adequate network code optimization and design. In fact. Proposition H] 
and Corollary |?] provide closed-form expressions of the coding gain for both scenarios where we account 
for realistic and ideal source-to-relay channels. A good criterion to design the network code, i.e., to exploit 
the inherent redundancy introduced by NC, might be to choose the generator matrix of the network code 
such that the following condition, for each source node, is satisfied: 

^0 Vi = 1,2,. . . , A^s and Vq = 1, 2, . . . , iVfl | (38) 

It is worth being mentioned that, in general, the most important criterion to satisfy is the diversity order 
requirement, as it has a more pronounced effect on the system performance. The optimization condition in 
( [38l ) can be taken into consideration if there is no reduction on the achievable diversity order for a given rate. 
Finally, we emphasize that both diversity order and coding gain can be adjusted by adding or removing relay 
nodes from the network, which, however, has an effect on the achievable rate as shown in Section |II1 The 
framework proposed in this paper can be exploited for many network optimizations, such as: i) designing 
the network code to achieve the best diversity order and coding gain for a given number of sources and 
relays {i.e., for a given rate); or ii) designing the network code to have the minimum number of source and 
relay nodes {i.e., to maximize the rate), for a given diversity order and coding gain. 

• EEP/UEP Capabilities. The diversity analysis in Proposition H] has pointed out that each source of the 
network can achieve a diversity order that is given by the separation vector of the network code. In other 
words, each source can achieve a different diversity order. In coding theory, this class of codes is known 
as UEP codes [49], and it can be very useful when different sources have to transmit data with a different 
quality-of-service requirement or priority. In other words, the network code might be designed to take into 
account the individual requirement of each source, instead of being designed by looking at the worst-case 
scenario only. For example, let us consider a network with three sources, with one of them having data to be 
transmitted with very low ABEP. Looking at the worst case scenario, we should optimize the system, and, 
thus, the network code as well, such that this source has, for a fixed transmit-power, a very high diversity 
order. If we cannot tune the diversity order of each source individually, we are forced to adopt a network 
code that provides the same high diversity order for all the sources of the network, which might have an 
impact on the achievable rate (see Section Our analysis showcases that UEP codes usually exploited in 
classical coding theory could be used to find the best trade-off between the diversity order achieved by each 
source and the rate of the network. In our opinion, this provides design flexibility, and introduces a finer 
level of granularity for system optimization, which has not been investigated yet for adaptive NC schemes. 
In fact, in general, network codes are designed such that all the sources have the same diversity order [23], 
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[27], [28]. Our framework provides a systematic way to guarantee unequal diversity orders for each source. 
Another interesting appUcation that might benefit from UEP capabilities provided by NC can be found in 
[65] and [66, Ch 3 and Ch. 5]. More specifically, in [65], UEP capability is called "incremental diversity". 
The idea is that energy consumption can be reduced if source nodes located farther from the destination 
can transmit with the same power as closer source nodes, and exploit UEP properties to achieve the same 
end-to-end performance. In other words, the incremental diversity offered by UEP network codes might 
be used to have even energy consumption among the nodes of the network with important implications for 
green applications [67]. Another application for energy saving is the exploitation of the proposed framework 
as an utility function for energy efficient network formation through coalition formation games [66, Ch. 5]. 

• Generalization of the performance analysis of dual— hop cooperative protocols. The framework proposed 
in this paper for can be thought as a generalization of the many results available in the literature for 
cooperative networks without NC. Among the many papers described in Section Jl let us consider, as an 
example, [16]. In [16], it is shown that a dual-hop three-node network using the DemF protocol can achieve 
full-diversity equal to 2 if the receiver has a reliable estimate of the instantaneous error probability at the 
relay. This result is included, as a byproduct, in our analysis, which is more general as it accounts for 
arbitrary sources, relays, and binary encoding vectors at each relay. In fact, under the classical coding theory 
framework, the distributed code used in [16] can be seen a repetition code with Hamming distance equal to 
2 for the single source of the network. Accordingly, from Proposition |?] we know that the diversity order is 
equal to 2, which confirms the analysis in [16] under a much broader perspective. In summary, the proposed 
framework can be used to study the end-to-end performance of dual-hop cooperative networks without NC, 
since a repetition code is a special network code. 

B. Impact of Receiver (Network) CSI on the Achievable Diversity 

In this section, we are interested in analyzing the importance of CSI at the receiver to achieve the full- 
diversity inherently available in the structure of the network code, which is given by its SV. In fact, it is 
important to emphasize that the conclusions drawn in Section [V-AI hold if the receiver has perfect knowledge 
of the cross-over probabilities computed in Section ITlI-AI This implies that the receiver knows the encoding 
vectors used at each relay node, along with the CSI of all the wireless links of the network. In general, the 
network code can be agreed during the initialization of the network or transmitted by each relay node over 
the control plane (at the cost of some overhead). On the other hand, CSI must be estimated at the receiver. 
In this section, we are aimed at showing the importance, to achieve full-diversity, of the knowledge of these 
cross-over probabilities. To this end, we assume that each receive node, including the destination, has access 
to the CSI of the wireless links that are directly connected to it (single-hop). In other words, the destination 
knows only the fading gains over the source-to-destination and relay-to-destination links, while it is not 
aware of the fading gains over the source-to-relay links. On the other hand, we assume that the destination 
is aware of the network code used at the relays. This is a requirement for any NC design. 
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With these assumptions, the destination is unable to compute the cross-over probabiUties in ([7]), and, thus, 
the received bits cannot be properly weighted according to their reliability. In such a worst-case scenario, 
the destination can only assign the same reliability to each received bit. This corresponds to set all the 
weights in Q equal to 1, i.e., w [m] = 1 for m = 1, 2, . . . , {Ns + Nf(). Accordingly, the demodulator in 
^ is no longer ML-optimum, and it simplifies to: 

fesi > fcSa . • • ■ : ''Sjv °c argmin {'^l^StD - bsA + ^ \bR,D - bnl > (39) 

i,Sie{o,i},-.-,6s„ge{o,i} [t=i ' 9=1' 'J 

By using the connection between network code design and classical coding theory described in Section 
|V-A[ the decoder in (|39l ) can be interpreted as a distributed Minimum Distance Decoder (MDD) applied to 
the overall network code [52]. The fundamental difference with classical coding theory is that, even though 
the receiver is not aware of CSI on the source-to-relay links, demodulation errors at the relay always 
take place and propagate through the network because of NC and forwarding operations. The demodulator 
in ( [39l ) simply cannot counteract these effects. Of course, this is a worst-case scenario as the destination 
has no estimates, even imperfect, of this CSI. The goal here is to understand the diversity order of this 
low-complexity but sub-optimal demodulator. Proposition \5\ provides an answer to this question. 

Proposition 5: Given the network model described in Section HH the demodulator in (l39l ) provides an 
end-to-end diversity order equal to (t = 1, 2, . . . , Ns): 



G'f''> = SV[t] ■ 



SV [t] 



(40) 



Proof: It follows by using the same steps as in Section|TV]by setting w [m] = 1 for m = 1, 2, . . . , {Ns + Nr). 
Due to space limitations, we describe only the main modifications of the proof that lead to (|40b . In particular, 
when w [m] = 1 for m = 1, 2, . . . , {Ns + Nr), and ^ simplify as follows: 

^ - n p E J n p [^[%,] ] (4^)"''"'^"'"' <5f^^ (41) 

keA [keA ) V 0/ 

(a) 

where: i) the large-SNR approximation in — )• is obtained by using the same development as in Appendix II-B I 



More specifically, in (1471) we have proved that Z can be seen as the error probability of a Maximum Ratio 
Combining (MRC) scheme with card {A} = dn (c, c) — d branches, where d = card {A}; and ii) G^J^'^ it 
related to the coding gain of T2, which is not shown here due to space limitations. By comparing (1451 ) and (|4TI ). 
it follows that T2 undergoes a reduction of the diversity order from du (c, c) to card {^} = du (c, c) — d. 
From (|27] ). because of the summation over d, each term of the APEP has no longer the same diversity order 
equal to du (c, c), but the allowed diversity orders fall in the range [du (c, c) — [dn (c, c)/2j , d// (c, c)]. 
Since end-to-end diversity is given by the addend having the smallest diversity order, we conclude that 
DIVapep = du (c, c) — max{d} = dn (c, c) — [du (c, c)/2j. Finally, by taking into account the relation 
between Hamming distance and SV given in Section |Vl (l40l ) is obtained. This concludes the proof. □ 
Proposition |5] brings to our attention the importance of the CSI of the source-to-relay links. In fact. 
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the demodulator in ( [39l ) loses approximately half of the potential diversity order inherently available in the 
network code. This result is in agreement with some studies available in the literature for simple cooperative 
networks without NC, such as [14], [16], and [22], where a similar diversity loss due to either non-coherent 
demodulation or imperfect CSI has been observed. Furthermore, this result seems to agree with the diversity 
that can be achieved by linear block codes over single-hop networks with hard-decision decoding [52, Sec. 
14—6-2]. However, it should be emphasized that in our case the diversity loss is not due to hard-decision 
demodulation at the physical layer, which is actually used for both demodulators in Q and ( [39l ). but it 
originates from the distributed nature of the network code, from demodulation errors at the relays, and from 
the demodulator that does not adapt itself to the reliability of the source-to-relay links. 

VI. Numerical and Simulation Results 

The aim of this section is to show some numerical examples to substantiate analytical derivations, claims, 
and conclusions of the paper. More specifically, we are interested in: i) showing the accuracy of the 
proposed framework for high-SNR, as well as the accuracy of diversity order and coding gain analysis; ii) 
understanding the impact of assuming ideal source-to-relay links, as it is often considered in the literature, 
and bringing to the attention of the reader that this might lead to misleading conclusions about the usefulness 
of NC over fading channels; iii) studying the impact of the network geometry on the end-to-end performance, 
and, more specifically, the role played by the positions of the relays; and iv) verifying the diversity reduction 
caused when the reliability of the source-to relay links is not properly taken into account at the destination. 
The analytical frameworks are compared to Monte Carlo simulations, which implements Q and ( |39l ) with 
no high-SNR approximations. Simulation parameters are summarized in the caption of each figure. 

a) Accuracy of the Framework for i.i.d. Fading Channels: Figs. [l]-[8] show the end-to-end ABEP for 
three network topologies (Ns = 2 and Nr = 2; Ns = 3 and Nji = 3; Ns = 2 and N^j = 5) and for different 
network codes. In particular, the network codes are chosen according to three criteria: i) NC is not used and 
only cooperation is exploited to improve the performance; ii) all the relay nodes implement binary NC on 
all the received data, as it is often assumed in the literature [25]; and iii) only some relay nodes perform 
NC on a subset of receiver packets. The first class of codes provides the reference scenario to understand 
the benefit of NC over classical cooperative protocols. The second class of codes represents the baseline 
scenario for network-coded cooperative networks. Finally, the third class of codes is important to highlight 
UEP capabilities, and to show that a non-negligible improvement can be obtained if the network code is 
properly designed and only some sources are network-coded. Numerical examples confirm the tightness of 
our framework for high-SNR, and that both diversity order and coding gain can be well estimated with 
our simple framework. Furthermore, the UEP behavior of many network codes can be observed as well. In 
particular, by comparing the SVs summarized in the caption of each figure with the slope of each curve, 
we can notice a perfect match, as predicted in Section |Vl Finally, we note that by comparing the results of 
the 2-source 2-relay network with the results of the 2-source 5-relay network, we can notice that if the 
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network code is not properly chosen, having multiple relays does not necessarily lead to a better diversity 
order. Since the rate of the system is smaller for larger networks (more relays), we can conclude that small 
networks with well-optimized network codes can outperform large networks where the network code is not 
adequately chosen. What really matters to optimize the performance of multi-source multi-relay networks 
is the SV of the network code, and, thus, the way the packets received at the relays are mixed together. 

b) Impact of the Source— to— Relay Links on the Achievable Performance: In Table IH we show a 
comparative study of the performance of three network topologies for realistic source-to-relay links, along 
with the scenario where cr|^^ — )• oo for t = 1, 2, . . . , Ns and g = 1, 2, . . . , Nr, which is denoted as "ideal" 
in the table. The results have been obtained from the analytical models and have been verified through Monte 
Carlo simulations. The accuracy between model and simulation for the "realistic" scenario can be verified 
in Figs. [l]-[8l since the same simulation setup is used. On the other hand, due to space limitations, similar 
curves for the "ideal" case are not shown, but similar accuracy has been obtained. The framework used for 
this latter scenario is given in Corollary |?1 As discussed in Section IIV-CI Table U confirms that there is no 
diversity loss between the two scenarios, but only a coding gain loss can be expected. This is because for 
both scenarios the ML-optimum demodulator is used. However, the conclusions about the usefulness of NC 
for both scenarios can be quite different. Let us consider, for example, the 2-source 2-relay network. In 
the "ideal" setting, there is no doubt that NC-3 and NC-4 should be preferred to NC-1 (no NC) and to 
NC-2 (all received data packets are network-coded), as one user achieves a higher diversity order while the 
other has the same ABEP as NC-1 and NC-2. On the other hand, the conclusion in the "realistic" setting 
is different. In this case, we observe that the higher diversity order achieved by one user is compensated by 
a coding gain loss for the second user. In other words, a coding/diversity gain tradeoff exists. However, this 
behavior is in the spirit of cooperative networking: one user might tolerate a performance degradation in a 
given communication round and wait for a reward during another communication round. Properly choosing 
the network code enables this possibility. Furthermore, by comparing NC-1 and NC-2, we can notice that 
different conclusions can be drawn about the usefulness of NC in the analyzed scenarios. In the "ideal" 
setting, a cooperative network with NC (NC-2) has the same ABEP as a cooperative network without NC 
(NC-1). The conclusion is that NC is useless in this case. On the other hand, the situation changes in the 
"realistic" setting. In this case, we can see that NC-2 is superior to NC-1, and, thus, we conclude that the 
redundancy introduced by NC can be efficiently exploited at the receiver when it operates in harsh fading 
scenarios. In fact, in the "realistic" setting, NC-2 can counteract the error propagation due to the dual-hop 
protocol, even though this network code is not strong enough to achieve a higher diversity order. Another 
contradictory behavior can be found when analyzing the 3-source 3-relay network. By comparing NC-1 
(no NC) and NC-2 (the relays apply NC to all received packets), we notice that in the "ideal" setting 
NC turns out to be harmful, as NC-2 provides worse performance than NC-1. On the other hand, in the 
"realistic" setting we notice that NC-1 and NC-2 provides the same ABEP. In other words, NC does not 
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help but at least it is not harmful. These examples, even though specific to particular networks and codes, 
clearly illustrate the importance of considering realistic source-to-relay links to draw sound conclusions 
about merits and demerits of NC for multi-source multi-relay networks over fading channels. Furthermore, 
we mention that, for all the network topologies studied in Table Jl NC-2 is representative of a network 
code that has been designed by keeping (l38l) in mind, as it provides the same high-SNR diversity order and 
coding gain for both "ideal" and "realistic" settings. Finally, we emphasize that our conclusions and trends 
depend on the coding gain of the network, whose study is often neglected due to its analytical intractability 
[16], [18], [19], [25]. In this paper, we succeeded to provide an accurate estimate of the coding gain as well. 

c) Accuracy of the Framework for i.n.i.d. Fading Channels and Impact of Relay Positions: In Fig. |9] 
and Fig. [TOl we analyze the accuracy of the framework for i.n.i.d. fading channels. We consider a 2-source 
2-relay network with nodes located as described in the caption of the figures. We consider five network 
topologies where the relay nodes can occupy different positions with respect to source and destination nodes. 
We observe a good accuracy of the framework, and notice that the positions of the relays can affect the 
end-to-end performance. This example shows that the proposed framework can be used, for arbitrary fading 
parameters, for performance optimization via optimal relay placement. 

d) Impact of Receiver CSI on the Diversity Order: In Fig. [TT] and Fig. [121 we study the impact 
of using the sub-optimal non-ML demodulator in ( |39l ). In particular, the ABEP of this demodulator is 
computed by using Monte Carlo simulations, and it is compared to the analytical investigation in Section 
IV-BI For comparison, the ABEP (analytical framework and Monte Carlo simulations) of the ML-optimum 
demodulator in Q is shown as well. The non-negligible drop of the diversity order can be observed, and, 
by direct inspection, it can be noticed that the curves have the slope predicted in (l40l ). This confirms the 
importance of CSI about the source-to-relay links in order to avoid substantial performance degradation. 

VII. Conclusion 

In this paper, we have proposed a new analytical framework to study the performance of multi-source 
multi-relay network-coded cooperative wireless networks for generic network topologies and binary encod- 
ing vectors. Our framework takes into account practical communication constraints, such as demodulation 
errors at the relay nodes and fading over all the wireless links. More specifically, closed-form expressions 
of the cross-over probability at each relay node are given, and end-to-end closed-form expressions of 
ABEP and diversity/coding gain are provided. Our analysis has pointed out that the achievable diversity 
of each source node coincides with the separation vector of the network code, which shows that NC can 
offer unequal diversity capabilities for different sources. Also, the importance of CSI about the source- 
to-relay channels has been studied, and it has been proved that half of the diversity might be lost if the 
reliability of the source-to-relay links is not properly taken into account at the destination. Monte Carlo 
simulations have been used to substantiate analytical modeling and theoretical findings for various network 
topologies and network codes. In particular, numerical examples have confirmed that the proposed framework 
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is asymptotically-tight for high SNRs. Finally, by comparing the performance of various network topologies, 
with and without taking into account decoding errors at the relays, we have shown that wrong conclusions 
about the effectiveness and potential gain of NC for cooperative networks might be drawn when network 
operations are oversimplified. This highlights the importance of studying the performance of network-coded 
cooperative wireless networks with practical communication constraints for a pragmatic assessment of the 
end-to-end performance and to enable the efficient optimization of these networks. The framework proposed 
in this paper provides an answer to this problem. 



Appendix I 

Proofs of Lemma[T1 LemmaIH and Lemma[3] 

A. Proof of Lemma |7] 

Lemma 1: Let Ti = Eh {nf=f '^^ P "^(c^c) 
and P [m], for m = 1, 2, . . . , Ns + Nr, given in Section Hill and in PropositionUl Then, over i.n.i.d. Rayleigh 
fading channels and for high-SNR, Ti has closed-form expression as follows: 



with P = [PsiD, • • • , -PSjvgD, PSi:NsRiD^ • • • ' PSi^nsPn^d] 
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where all symbols are defined in Proposition \3\ 

Proof: Owing to the assumption of independent fading channels, it follows, by direct inspection, that 
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H for m = 1,2, . . . , + iV/j are independent RVs, and, thus, Ti = nfc=i'''''^ Eh |p m^^%^ | = 

Eh |p "i-l^^-j |. Furthermore, from the definition of P [m] in 
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Section JII] and Proposition |7] for high-SNR we have 
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The results in (1431) can be obtained from the following chain of equalities and high-SNR approximations: 
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where: i) = comes from [52, Eq. (14-3-7)]; ii) — >■ is the high-SNR approximation of = in [52, Eq. 

(62) (fei) 
(14-3-13)]; iii) — )• is the high-SNR approximation of = , which simply neglects the term Psi.^ r^Pr^d, 

(ci) 

as it decays faster for high-SNR; iv) — >• follows by noticing that 1 — ^gs^R^Ps^Rg — ^ 1 for high-SNR; and 

(C2) (02) (ci) . . . . 

v) — )• is, similar to — > , is the high-SNR approximation of — >• . From (1441 ). (142b follows by using notation 

and vector representation in Proposition 13 More specifically, the vector Ac,c takes into account that only the 

— — — (G) 

indexes in the set © (c, c) = {m\c [m] / c [m]} have to be included in Ti, and the vector Sgj^j-, accounts 

for the dual-hop relaying protocol and the specific network code. This concludes the proof. □ 
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Two important remarks are worth being made about Lemma [7] First, we would like to emphasize that, 
for ease of presentation and to stay focused on the most important issues of our analysis, i.e., dual-hop 
networking and NC, the results in (l42l) and (|43] ) are here given for Rayleigh fading only. However, they can 
be generalized to other fading distributions for which the high-SNR approximation in [61] exists. In this 
paper, Rayleigh fading is studied for illustrative purposes only. Second, by comparing — )• and [62, Eq. (40)], 
it follows that, for high-SNR, the effect, on the error probability at the relays, of performing NC on noisy 
and faded received data is equivalent to an Amplify-and-Forward (AF) relay protocol with CSI-assisted 
relaying [60] and with a number of hops equal to the number of sources that are network-coded at each relay. 
This conclusion is in agreement with the equivalence between the error probability at the relays and the error 
performance of DemF relay protocols already highlighted in Section ITlI-AI In fact, in [58] it has been shown 
that, except when the number of hops is very large and the fading severity is very small, the performance 
of AF and DemF protocols is very close, for high-SNR, to each other. As the number of sources that can 
be network-coded is, for practical applications, not very large, this high-SNR approximation can be very 
useful to get formulas that provide insights on the system behavior. The high-SNR equivalency between 
( |43l ) and AF relaying is exploited in Lemma |2] to get high-SNR but closed-form and accurate formulas. 



B. Proof of Lemma |2] 

Lemma 2: Let us consider the term T2 = Eh |min |]^ 
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and in PropositionUl Then, over i.n.i.d. Rayleigh fading and for high-SNR, T2 has closed-form expression: 
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where all symbols are defined in Proposition \3\ 

Proof: The computation of T2 is very analytically involving. To get accurate, but closed-form and insightful 
formulas that can shed lights on the network behavior, we exploit some high-SNR approximations. More 
specifically, the starting point is the following high-SNR approximation: 

PStD = Q (^^2 {E^/No)\hs,D\^^ 

(a) I 
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The approximation in — )■ follows from the closing comment in Lemma |7] where we have proved that 
for high-SNR the cumulative error due to performing NC on wrong demodulated bits at the relays can 
be well-approximated by an equivalent AF multi-hop relay network with a number of hops that is equal 
to the number of network-coded sources. In particular, from [60] we can recognize that the argument of 

(a) 

the Q-function in — ^ is the end-to-end SNR of an AF relay network, which takes into account the relay- 
to-destination link and the cumulative error due to combining, at the most, Ns source. The number of 
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sources that are actually network-coded depends on the number of non-zero NC coefficients gstRg- Let 
us emphasize that the formula for the direct source-to-destination link is exact, but we have decided to 
re-write it to better understand that high-SNR approximation applies only to the signals forwarded from 
the relays. Thus, in T2, we have P ^ Q (^2 (S„/iVo) SNRj^^y , where: i) SNr|^^)_) = {hs^nf 

for the source-to-destination Unks, and by bearing in mind that in this case we have a true equality; and ii) 

for the relay-to-destination Unks. Thus, T2 simplifies: 
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where: i) the approximation in —5- is proved in Lemma \3\ ii) Ti and T2 are two constant factors whose 

1-1 (^) 
closed-form expression is given in Lemma |2J and iii) the equality in = comes from the fact that the 

Q-function is monotonically decreasing for increasing values of its argument. 

The last expression in (|47] ) has a convenient structure that can be averaged over fading channel statistics. 

To this end, the following considerations can be made: i) T2 can be seen as the ABEP of a dual-branch 

Selection Combining (SC) scheme, where the equivalent SNR of first and second branch is [53] SNRi = 

2(^„/iVo)TiEfee^SNR|^J_) and SNR2 = 2 (S„/iVo) T2 Efce^ SNr|^)_), respectively; ii) both SNRi 

and SNR2 can be seen as the equivalent SNR of a Maximum Ratio Combining (MRC) scheme with a 

number of branches given by card{^} = d and card{^} = dn {c,c) — d, respectively; and iii) the 

"virtual" SC and MRC branches contain independent RVs, as it can be verified via direct inspection. Thus, 

a closed-form and high-SNR approximation of T2 in (|47] ) can be obtained by using the method in [61]. 

More specifically, by considering: i) the definition of SNr|^ in ( |47l ); ii) the closed-form expressions of 

Ti and T2 in Lemma \3l and iii) the general parametrization in [61, Prop. 1, Prop. 4] for systems with 

receive-diversity, we can obtain, after lengthly algebraic manipulations, the final result shown in (1451 ). In 

particular, Ac,c and Sg^j-, have the same meaning as in Lemma [7J while the term into the square brackets 

(b) 

accounts for the SC/MRC high-SNR approximation of = in (1471 ). This concludes the proof. □ 
Finally, similar to Lemma [7] we emphasize once again that the closed-form solution in (1451 ) can be 
generalized to other fading channel models by using [61] and [62]. 
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C. Proof of Lemma\3\ 

Lemma 3: Let A = Ilfce^ Q {Em/ No) SNRj^^ with SNRj^^_^ defined in Lemma ^ Then, for 
high-SNR and Rayleigh fading, A can be tightly approximated as follows: 



^ = n '3 ( V ^ (i?m/A^o) SNR|;:;_j J ^ Q /2 (E^/No) T ^ SNR 
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and d = card {A}. 



where T — 

Proof: From the Chemoff bound, i.e., Q (x) < (1/2) exp (— a;^/2) < exp (— which is accurate for 
X :$> 1 that in our case implies high-SNR (Em/No » 1), the following approximation holds: 
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where T is a constant correction term, which is introduced to recover the coding gain inaccuracy that might 
arise when using the Chernoff bound [16]. The high-SNR approximation in (l49l ) can be explained as follows. 
By direct inspection, left- and right-hand side terms can be shown to have both diversity order equal to 
d = card{^}. In fact, the left-hand side is the product of d terms each one having diversity one. On the 
other hand, the right-hand side term is the error probability of a MRC scheme [53] with d diversity branches 
at the receiver, which is known to have diversity d [61]. The constant (correction) factor T is introduced 
only to avoid coding gain inaccuracies, which are always present when using the Chernoff bound. Since the 
goal of this paper is to accurately estimate both coding gain and diversity order, the accurate evaluation of 
T is instrumental to estimate the end-to-end performance of the system. 

To get an accurate, but simple and useful for further analysis, approximation we use first-order moment 
matching to estimate T in (|49l ). The motivation is that, as we will better substantiate at the end of this proof, 
it allows us to have a closed-form estimate of T that depends only on d in (1481 ). while it is independent of 
the fading parameters. In formulas, we seek to find T such that the following equality is satisfied: 
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To this end, we need closed-form expressions of both averages in ([50] ). Once again, we use the high-SNR 
parametrization in [61], which leads to the following result: 
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where: 1) — )• is obtained by taking into account that (i) SNR; are statistically independent for k & A; 
(ii) according to Lemma |2J SNR^ , can be seen as the end-to-end SNRs of an equivalent multi-hop AF 
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relay protocol; and (iri) by using asymptotic analysis for multi-hop AF relay networks in [62]; and 2) 
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is obtained from [61] and [62] by recognizing that we have to compute the average of an equivalent MRC 
scheme where each branch is an equivalent multi-hop network that uses the AF relay protocol. Finally, by 
equating the two terms in dST] ). T in ( |48] ) can be obtained. As mentioned above, T is independent of channel 
statistics. Similar to Lemma [7] and Lemma |2] we mention that the proposed procedure can be applied to any 
fading channel model, for which the parametrization in [61] is available. This concludes the proof. □ 
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ABEPgj = yGi^*\rnj . WHERE G^f*^ IS THE DIVERSITY ORDER, Gc'^*^ IS THE CODING GAIN, AND WE HAVE DEFINED 

7m = Bm/A^O- 2-SOURCE,2-RELAY NETWORK. NC-1: fefl^ =6Sifli,f'fl2 = 6S2R2 ^ NC-2: &Rj =f>SiBi e^SsHi. 

bR2 =bSifl2 ©^S2fl2;NC-3:6fl^ = i>SiRi ®is-,Ri,bR-, = fcs2fl2 ^ NC-4: 6^^ =6sifli.bi{2 = 6sifl2 ® ''S2i?2 • 3-SOURCE, 3-RELAY 
NETWORK. NC- 1 : bfli = , bii2 = H2 , ''Ha = fl3 ; NC-2 : &Ri = e Sg^ jii ffi 653 jj^ , bjjj = Sgj H2 © &S2 H2 © ^Sg fla - 

6j{3 =6SiH3©&S2J{3ffi6s3H3;NC-3:6Hi =Ss^Hi.6i{2 =&S2H2'''fl3 = bsiHg ©6s2fl3 ffiSs3il3-2-SOURCE,5-RELAY NETWORK. 
NC-l:6flj = i>SiRi,bR2 =6siR2'*fl3 =''SiH3.fefl4 =''S2H4.^H4 = 6s2Hr, ' '^'^"2: feflj = bSifli © 632^1 > ^^2 = ^SiH2 ® *S2 fl2 > 

bj{3 = &Sifl3 ffi i>S2R3, *fl4 = ''Sifl4 © &S2fl4. ''fls = ''Sifls © i>S2R5< NC-3: fcflj^ = bSifli. i>R2 = 6Sifl2' ''fl3 = bSiRs ® bS2R3' 

bRi = bsiRi © &S2H4 . 6j{5 = bs2R5 ■ Finally, i.i.d. fading with o-g = 1 is considered. 
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Network: 3-source, 3-relay 
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Network: 2-source, 5-relay 
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Fig. 1. ABEP of a 2-source 2-relay network. Markers show Monte 
Carlo simulations and solid lines show the analytical framework. Setup: 
i) i.i.d. fading with erg = 1; and ii) bji^ = bsiRi ffi ''Safli. ''^2 = 
^SiR2 ® bs2R2 - Separation Vector is SV = [2, 2]. 
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Fig. 2. ABEP of a 2-source 2-relay network. Markers show Monte 
Carlo simulations and solid lines show the analytical framework. Setup: i) 
i.i.d. fading with erg = 1; andii) bn^ = bsiRi,bii^ = fesi^a ®^S2ii2 ■ 
The Separation Vector is SV = [3, 2]. 
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Fig. 3. ABEP of a 3-source 3-relay network. Markers show Monte 
Carlo simulations and solid lines show the analytical framework. Setup: 
i) i.i.d. fading with erg = 1; and ii) bn^ = bs^Rj^, bji^ = bs2fl2' 
bflg = . The Separation Vector is SV = [2, 2, 2]. 
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Fig. 4. ABEP of a 3-source 3-relay network. Markers show Monte 
Carlo simulations and solid lines show the analytical framework. Setup: 
i) i.i.d. fading with erg = 1; and ii) bR^ = bSiRi ffi ^S2iii © ^SsRi, 

bR2 = ^SiR2 ffi ^S2fl2 © ^S3fl2' ^flj = ^Siflj © ''52^3 © ^SsRa- 

The Separation Vector is SV = [2, 2, 2]. 
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Fig. 5. ABEP of a 3-source 3-relay network. Markers show Monte 
Carlo simulations and solid lines show the analytical framework. Setup: i) 
i.i.d. fading with erg = 1; and ii) b^^ = SsiKj, bn^ = 632^2' ^^3 = 
^SiRs © ''S2-R3 © ^SsR^- The Separation Vector is SV = [3, 3, 2]. 
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Fig. 6. ABEP of a 2-source 5-relay network. Markers show Monte 
Carlo simulations and solid lines show the analytical framework. Setup: 
i) i.i.d. fading with erg = 1; a.nd ii) b^^ = bsj^Rj^, bu^ = 6siiJ2' 
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Fig. 7. ABEP of a 2-source 5-relay network. Markers show Monte 
Carlo simulations and solid lines show the analytical framework. Setup: 
i) i.i.d. fading with erg = 1; and ii) ^jji = ^Sifli © ''Safli' = 
bSiR2 ffi^S2H2'.^«3 = ''S1H3 ® bs^Rs, = bs^R^ ® bs^R^, 

bijj = feSi-Rs ® bs2R5 - Separation Vector is SV = [2, 2]. 
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Fig. 8. ABEP of a 2-source 5-relay network. Markers show Monte 
Carlo simulations and solid lines show the analytical framework. Setup: 
i) i.i.d. fading with erg = 1; and ii) bjii = bsifli, ''fl2,= ^Sifl2' 
bRs = bsiRs ® ''S2H3. bft^ = f)Sifl4 ffi f'S2-R4. ^-Rs = bs^R^- The 
Separation Vector is SV = [5, 4]. 
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Fig. 9. ABEP of a 2-source 2-relay network. Markers show Monte 
Carlo simulations and solid lines show the analytical framework. Setup, 
i) i.n.i.d. scenario with a = 3, ct\y ~ '^XY' ''^^ nodes are located 
at positions (in meters): Si = (0,25), ^2 = (0,-25), D = (50,0), 
Ri = (xjj^,12.5), R2 = (xh2,-12.5); and iii) b^^ = bg-^R^ ® 
bs2Ri' bR2 = bs2R2 - Furthermore, we have: i) xji^ = 25 and x^^ = 
25 in Scenario I; ii) xn-^^ = 5 and xn^ = 5 in Scenaiio 2; iii) xn-^ = 45 
and xr„ = 45 in Scenario 3; iv) xr, = 5 and xn„ =45 in Scenario 
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Fig. 10. ABEP of a 2-source 2-relay network. Mai'kers show Monte 
Carlo simulations and solid lines show the analytical framework. Setup, 
i) i.n.i.d. scenario with a = 3, cr^y = '^XY^ ''^ "'^ nodes are located 
at positions (in meters): Si = (0,25), S2 = (0,-25), D = (50,0), 
^1 = [xji^, 12.5), R2 = (xfl2,-12.5); and iii) b^^ = Ssiflj. 
bji^ = bs2R2 - Furthermore, we have: i) xn-^ = 25 and xji^ = 25 in 
Scenario 1; ii) xn^ = 5 and xn^ = 5 in Scenario 2; iii) xn-^ = 45 
and = 45 in Scenario 3; iv) xji-^ = 5 and x^, = 45 in Scenario 
4; v) xn-^ = 45 and xn^ = 5 in Scenario 5. 
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Fig. 1 1 . ABEP of a 3-source 3-relay network. Markers show 
Monte Carlo simulations for the ML demodulator, dashed lines show 
the analytical framework for the ML demodulator, and solid lines with 
mai'kers show Monte Carlo simulation of the non-ML demodulator. 
Setup: i) i.i.d. fading with o-^ = l;andii)6flj = bsiiii.fefl2 = ''S2fl2' 
fcfla = '^SiRs © bs2R3 ® bSsRs- 
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Fig. 12. ABEP of a 2-source 5-relay network. Markers show 
Monte Carlo simulations for the ML demodulator, dashed lines show 
the analytical framework for the ML demodulator, and solid lines with 
mai'kers show Monte Carlo simulation of the non-ML demodulator. 
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